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Abstract. A microscopic magnetic model for the spin-1/2 Heisenberg chain 
compound CuSe205 is developed based on the results of a joint experimental and 
theoretical study. Magnetic susceptibility and specific heat data give evidence for quasi- 
1D magnetism with leading antiferromagnetic (AFM) couplings and an AFM ordering 
temperature of 17 K. For microscopic insight, full-potential DFT calculations within 
the local density approximation (LDA) were performed. Using the resulting band 
structure, a consistent set of transfer integrals for an effective one-band tight-binding 
model was obtained. Electronic correlations were treated on a mean-field level starting 
from LDA (LSDA+J7 method) and on a model level (Hubbard model). In excellent 
agreement of experiment and theory, we find that only two couplings in CuSe205 are 
relevant: the nearest-neighbour intra-chain interaction of 165 K and a non-frustrated 
inter-chain coupling of 20 K. From a comparison with structurally related systems 
(Sr2Cu(P04)2, Bi2Cu04), general implications for a magnetic ordering in presence of 
inter-chain frustration are made. 
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1. Introduction 

Low-dimensional spin-1/2 systems attract much interest due to a variety of ground 
states (GS) found in these systems which originates from an interplay between different 
exchange interactions and strong quantum fluctuations. There are, for instance, the 
spin-Peierls GS in CuGe03 pL|, the helical GS in LiCuV04 [2], and the quantum 
critical behaviour in Li 2 ZrCu04 [3j H] etc. Besides, many of these materials (mostly 
cuprates, vanadates and titanates) appeared to be realizations of theoretically long- 
studied models in good approximation. One of the most prominent models is the spin- 
1/2 nearest-neighbour (NN) chain described by the Heisenberg Hamiltonian, for which 
the exact solution has been derived by Bethe [5]. The first compounds proposed to be 
good material realizations of this model were Sr 2 Cu0 3 and Ca 2 Cu0 3 [6j [7J. Recently, 
Sr 2 Cu(P04)2 and Ba 2 Cu(P04) 2 were suggested as even better realizations [8] [12]. 
followed by a study of K 2 CuP 2 7 [13] that qualified this compound to be the best 
realization of the spin-1/2 NN Heisenberg chain up to date. As a natural consequence 
of its simplicity, this model poorly describes one-dimensional and quasi-one-dimensional 
systems where additional interactions, like longer range couplings or anisotropies, are 
present. Thus, extensions of this model are required to allow an accurate description of 
real materials. 

The simplest extension of the model is the inclusion of a next-nearest-neighbour 
(NNN) coupling J 2 leading to the so-called zig-zag chain model. In case of an 
antiferromagnetic (AFM) J 2 both NN and NNN couplings cannot be simultaneously 
satisfied, in other words, the system is magnetically frustrated. Here, the intra-chain 
frustration enriches the phase diagram with the spiral GS, the gapped AFM GS and a 
quantum critical point at J 2 /Ji = —0.25 [2] [UJ. The evaluation of the two parameters 
in the zig-zag chain model allows to estimate quantities which can be directly measured 
or derived from experiments, namely spin-spin correlation functions, thermodynamic 
properties and the response in high magnetic fields. Nevertheless, this model fails to 
describe phenomena like long-range magnetic ordering, since one- or two-dimensional 
systems do not order at finite temperatures according to the Mermin- Wagner theorem. 
Thus, to account for magnetic ordering, the inter-chain (IC) coupling has to be included 
in the model. This problem has been addressed in a series of theoretical works [18]- 
[25] . but the simplifications that had to be made to keep the models solvable (at 
least approximately) impede an accurate description of complex situations. The spin- 
1/2 Heisenberg chain system Sr 2 Cu(P04) 2 is the most prominent example for a huge 
discrepancy (two orders of magnitude) between the theoretical (adopting a simple NN 
inter-chain coupling) and the experimentally observed ordering temperature. The origin 
of this discrepancy is hidden in the effects of magnetic frustration and anisotropy, and 
the disentanglement of the effects is difficult. Thus, a reliable theoretical description of 
this phenomenon is still missing. 

A natural way towards a deeper understanding is the search for real material 
realizations of "easy" models. For such systems, the experimental data can be 
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supplemented by reliable microscopic models. This way, joint experimental and 
theoretical studies can challenge and improve the existing theoretical approaches. 

Cu 2+ phosphates are up to date the best realizations of the NN spin- 1/2 Heisenberg 
chain model [U [TOj [13] and could have a great potential for the discovery of further 
low-dimensional systems. Unfortunately, the experimental information about these 
materials is rather limited since they are up to now available as powders only, 
although several attempts have been made to grow single crystals required for advanced 
experimental studies. 

In this paper, the structurally closely related compound, CuSe 2 5 , is investigated. 
Since selenites are often susceptible to chemical transport, the advantage of this material 
is the potential to grow large single crystals of high quality. Previous studies on a powder 
sample [26] may hint a ID character of its magnetic properties, but the low-temperature 
data are strongly affected by impurities (figure 3 in |26j). Therefore, to probe the ID 
nature of the system, a new detailed study on high quality samples with lower defect 
concentration is desirable. 

Though the chain-like arrangement of CuC>4 squares in CuSe 2 C>5 is topologically 
similar to that in Cu 2+ phosphates, the geometry of magnetic coupling paths between 
the structural chains is essentially different. Thus, the role of magnetic frustration, 
which is ruled by the inter-chain coupling, can be evaluated in a comparative study. 

The paper is organized as follows. In section [2] we describe the synthesis, sample 
characterization and experimental as well as theoretical methods used in this work. 
In section EJ we discuss the crystal structure of CuSe 2 Os in comparison to related 
systems. Section @] reports the results of our measurements and theoretical calculations 
and proposes an appropriate microscopic model. A brief summary and an outlook are 
given in section [5j 

2. Method and sample characterization 

Single crystals of CuSe 2 5 were grown by chemical vapour transport using TeCl 4 as a 
transport agent. Using a micro-crystalline powder of CuSe 2 Os (obtained from a mixture 
of CuO and Se0 2 at 723 K) as a source, the transport experiments were carried out in 
an endothermic reaction of T 2 (source) 653 K to T x (sink) 553 K. 

The obtained crystals have a green colour and form strongly elongated (along [001]) 
plates, which macroscopically look like needles. The typical length of a needle is 5- 
10 mm and the width does not exceed 1 mm and for most crystallites it is considerably 
smaller. The slight disorientation of plates forming a needle intricates a precise X-ray 
diffraction measurement on single crystals. Thus, the samples were characterized by X- 
ray powder diffraction and energy- dispersive X-ray spectroscopy (EDXS) experiments. 
The lattice parameters of the synthesized crystals are similar to those reported for 
CuSe 2 5 (tabled]). The results of the EDXS analysis (Cu 32.78±0.31, Se 67.14±0.23) 
for 13 points (2 crystals) yield Cu:Se ~ 0.488±0.006, very close to the ideal ratio of 0.5. 
Thus, the obtained single crystals represent an almost pure CuSe 2 Os phase. 
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Magnetization was measured in a SQUID magnetometer (1.8 - 350 K) in magnetic 
fields up to 1 T. Heat capacity (1.8 - 100 K) was determined by a relaxation method 
up to hqH = 9 T. 

DFT calculations were carried out using the full potential local orbital code (FPLO) 
version 7.00-27 [27 \. The standard basis set and the Perdew-Wang parameterization of 
the exchange-correlation potential were used (28]. Strong on-site Coulomb interaction 
in the Cu 3d orbitals, insufficiently described in the LDA, was taken into account 
independently (i) by mapping the LDA antibonding Cu — O dp a bands onto a tight- 
binding model (H = J2i e i + J2<ij>a Uj( c icr c j,a + H.c.)) and subsequently via a Hubbard 
model (H = J2i e i + J2<ij>a Uj(ci,aCj jrT + H.c.) + U e s J2i "<i,T ^i>l) onto a Heisenberg model 
(H = J2<ij> JijSiSj) (the procedure is well justified for spin excitations in the strongly 
correlated limit (Ucs ^> Uj) at half-filling ((rij) = 1)) and (ii) by using the LSDA+t/ 
method [29] (Ud = 6.5 eV, Jd — 1 eV). For the LDA calculations, we used a /c-mesh 
of 1296 Appoints (355 points in the irreducible wedge), for LSDA+£7 calculations of 
supercells irreducible fc-meshes of 226, 242, 147 and 126 /c-points were used. All k- 
meshes are well converged. 

Quantum Monte-Carlo simulations have been preformed on N = 1200 sites clusters 
of S = 1/2 spins (30 coupled chains of 40 sites each) using the ALPS software 
package [30J. 

3. Crystal structure and empirical magnetic models 

Crystal structures of cuprates are often subdivided into four large groups according 
to their dimensionality, which reflects how their elementary building blocks — CUO4 
plaquettes (planar or distorted) are connected: they can be isolated (zero- 
dimensional, 0D) or form chains (one-dimensional, ID), layers (two-dimensional, 2D) 
or frameworks (three-dimensional, 3D). Although it is true for many systems that 
the magnetic dimensionality follows the dimensionality of the crystal structure, the 
real situations are often more complex, especially for 0D cases. There, the magnetic 
dimensionality is ruled by (i) the orientation of neighbouring plaquettes, and (ii) the 
position of anion groups formed by non-magnetic atoms that bridge the magnetic 
plaquettes. In most cases, the connection between structural peculiarities and the 
appropriate magnetic model cannot be accounted for by applying simple empirical rules 
(for instance, Goodenough-Kanamori-Anderson rules). 

Therefore, an almost complete understanding of the macroscopic magnetic 
behaviour for a certain system of this class can be achieved only basing on a 
relevant microscopic model. The latter can be constructed either by using advanced 
experimental techniques (for instance, inelastic neutron scattering) or theoretical (DFT) 
calculations. Naturally, the most reliable approach is the combination of such a theory 
and experiment. Due to the complexity of such an analysis, it has been accomplished 
only for a rather limited number of real systems. 

Two well studied systems of this class - - Bi 2 Cu04 [31] and Sr 2 Cu(P04)2 [8] 
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- are both structurally OD cuprates, but antipodes with respect to their magnetic 
behaviour. A drastic change of the magnetic coupling regime originates from the 
arrangement of neighbouring plaquettes (figure [Q, right panel: top and bottom): 
stacking (accompanied by additional twisting) of neighbouring plaquettes on top of 
each other makes Bi 2 Cu0 4 a 3D magnet with Tn ~ 47 K [32J while in Sr2Cu(P04)2 
the plaquettes are arranged in a planar fashion (formally reminiscent of an edge-sharing 
chain with every second plaquette cut out), leading to a pronounced ID behaviour and 
a very low Neel temperature Tn = 0.085 K [9]. In this context we mention that the 
magnetic dimensionality of these systems is controlled by the dihedral angle <fi between 
neighbouring plaquettes (figure CD, right panel, middle). In CuSe205, the magnetic 
plaquettes are isolated (like in Bi 2 Cu0 4 and Sr 2 Cu(P0 4 )2) but tilted with respect to 
each other forming a dihedral angle <fi of about 64°, i.e. in between = 0° for the 
3D Bi 2 Cu0 4 and = 180° for the ID Sr 2 Cu(P0 4 ) 2 (figure GQ right panel). Another 
controlling parameter is the direct Cu — Cu distance d. Again, CuSe 2 5 with d close to 4 
A lies in between d = 2.9 A for Bi 2 Cu0 4 and d — 5.1 A for Sr 2 Cu(P0 4 ) 2 . Thus, CuSe 2 05 
is structurally in between the two closely related systems: the 3D magnet Bi 2 Cu0 4 and 
the ID magnet SrCu 2 (P0 4 ) 2 . Does this analogy hold also for the magnetism? To answer 
this question, additional arguments have to be addressed. 

Besides and d, further structural features provide a deeper insight into the crystal 
chemical aspects relevant for the magnetism. In CuSe 2 Os, two SeC>3 pyramids sharing 
an oxygen atom (forming Se 2 5 polyanions) bridge neighbouring Cu0 4 plaquettes (see 
figure [I]) . This structural peculiarity is reflected in the morphology of the synthesized 
crystals (section [2]): the needle-like shape with an elongation along [001] fits perfectly to 
structural chains along c (figured]) formed by alternation of Cu0 4 plaquettes and Se 2 Os 
polyanion groups. From the topological similarity of this structural chain to the one in 
the structure of Sr 2 Cu(P0 4 ) 2 (there, neighbouring plaquettes are bridged by two P0 4 
tetrahedra) a ID behaviour of CuSe 2 Os might be expected. A second argument supports 
this proposition: in CuSe 2 Os, the structural chains are not connected by covalent 
bonds, making a strong inter-chain coupling unlikely, similar to Sr 2 Cu(P0 4 ) 2 , where 
the neighbouring chains are well separated by Sr cations. These similarities of CuSe 2 Os 
and Sr 2 Cu(P0 4 ) 2 may lead to the conclusion that both systems imply essentially the 
same physics. However, a closer inspection of more subtle crystal chemical aspects 
immediately reveals an important difference related to the inter-chain coupling. As 
it follows from the microscopic model [8], in Sr 2 Cu(P0 4 ) 2 there are two relevant NN 
inter-chain couplings (2.7 K), which are equivalent by symmetry. Together with an intra- 
chain NN coupling (187 K), they induce magnetic frustration which commonly leads to a 
considerable decrease of the ordering temperature (T N =0.085 K). In CuSe 2 5 , these two 
inter-chain couplings are not symmetry-equivalent. Therefore, by reducing the strength 
of one of them, the frustration can be lifted. 

So far, crystal chemical considerations provided us with a qualitative insight. For 
a quantitative model, a microscopic analysis is required. Thus, in the next section a 
microscopic model basing on the results of DFT calculations is constructed. 
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Figure 1. Left panel: the crystal structure of CuSc2 0s. Isolated CuC>4 plaquettcs 
(yellow) are bridged by SeC>3 pyramids (gray) and form chains running along c. The 
chains are closely stacked in b direction, and well separated in the a direction. Right 
panel: geometry of "chains" formed by isolated CUO4 plaquettes. Neighbouring 
plaquettes are stacked and twisted with respect to each other in Bi 2 Cu04 (top), tilted 
in CuSe2 05 (middle) and form planar edge-sharing chains with every second plaquette 
cut out in Sr2Cu(PC>4)2 (bottom). The non-magnetic groups (BiC>4, SeC>3 and PO4 
for Bi2Cu04, CuSe20.5 and Sr2Cu(PC>4)2, respectively) bridging the neighbouring 
plaquettes are not shown. 



A prerequisite for an accurate modeling based on a band structure code is reliable 
structural information. For CuSe205, two refinements of the same structural model 
(space group C2/c with four formula units per cell) have been proposed so far [331 EH]- 
Both structural data sets agree quite well with each other and with the lattice parameters 
of the synthesized samples (see table [Tj) Q The reliability of the structural data has been 
indirectly confirmed a posteriori by the good agreement of calculated and experimentally 
measured quantities. 

4. Results and discussion 

4-1. Thermodynamical measurements 

The first probe for the magnetic properties of a certain system is the measurement 
of magnetization (M) at various temperatures in a constant field (H) yielding 

f: It is well known, that X-ray diffraction analyses may result in considerable inaccuracies for internal 
coordinates of light elements (especially, hydrogen). These inaccuracies can have a large impact on the 
magnetic properties |35j . Since there are no light atoms in CuSe205, we rely on the diffraction analysis 
and therefore, no structural relaxation has been performed. 
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Table 1. Comparison of measured lattice parameters a, b, c, the monoclinic angle /3 
and the unit cell volume V of CuSe20s with previously published data. 



parameter Ref. [33] Ref. [33J this work 



a, A 


12.3869 


12.254 


12.272 


b, A 


4.8699 


4.858 


4.856 


c, A 


7.9917 


7.960 


7.975 


/?, ° 


109.53 


110.70 


110.91 


v, A 3 


447.13 


443.27 


443.95 



the temperature dependence of magnetic susceptibility (x(T) = M{T)/H). This 
measurement already yields valuable information on the magnetic dimensionality, the 
sign and the energy scale of leading couplings, the presence of a spin gap, the spin 
anisotropy and the quality (defects, purity) of a sample. 




T(K) T(K) 



Figure 2. Left panel: magnetic susceptibility of CUSC2O5 as a function of 
temperature. The magnetizing field is 10 kOe. For the graphic presentation, we show 
only one of each five measured points. The Bethe ansatz fits arc shown with dashed and 
solid lines. Right panel: inverse magnetic susceptibility as a function of temperature. 
Curie- Weiss fits (for T > 230 K) are shown with lines. The temperature-independent 
contribution \o i n the Curie- Weiss fits was set to zero. 

For CuSe 2 5 , the magnetic susceptibility curves for both field orientations (figure [2j 
left panel) have a broad maximum at T max m 101 K and a finite value of x a t the 
lowest temperature measured (1.8 K), indicating the low-dimensional behaviour and 
the absence of a spin gap. The high-temperature parts of the curves obey the Curie- 
Weiss law X (T) = C/(T + 8) (figure^ right panel; T > 220 K, Hf 9 = 165 K, C = 0.51, 
g = 2.32; H ± : 9 = 170 K, C = 0.43, g = 2.15). The positive Curie- Weiss temperature 
evidences that the dominating couplings in CuSe20s are antiferromagnetic. The shape 
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of the experimental curve reveals a close similarity to a spin-1/2 Heisenberg chain 
model. This model has an exact solution given by Bethe ansatz [5] and parametrized 
by Johnston et al. [36]. We have fitted the experimental curves using the parametrized 
solution (figure El left panel). 

To account for the deviation of the fitted curves from experimental ones we have 



the magnetic susceptibility measured perpendicular to needle-like crystallites can be 
perfectly fitted by a consistent set of parameters (J = 157.6 K, g = 2.00) in the whole 
temperature range down to the ordering temperature, while the fit to the susceptibility 
measured parallel to the chains (J = 157.1 K, g — 2.14) shows deviations at low 
temperatures (below T max ). This difference likely originates from a slight misalignment 
of microscopic plates in the needle-like crystallites. 

A phase transition is observed at 17 K for both orientations of the magnetizing 
field. The nature of this magnetic transition can be understood by examination of 
the low-temperature part of the curve (below the transition). The interpretation 
is straightforward as soon as we account for (i) impurity effects and (ii) effects of 
misalignment of the sample (relevant especially for H\\, as shown above). Due to the 
high quality of samples, the temperature region between the kink at 17 K down to 10 K 
is practically unaffected by defects (no Curie tail). In this range, x± decreases very 
slightly, while x\\ drops distinctly on cooling, following the theoretical result for ordered 
collinear antiferromagnets [37J. We assign the slight decrease of x± (theory predicts it 
to be constant) to a small misalignment of crystallites, in agreement with the deviations 
of the Bethe ansatz fit. The small upturn in x\\ & t about 6 K (a zero susceptibility 
at zero temperature follows from theory) is likely related to defects and paramagnetic 
impurities. 

To get additional information about the magnetic properties, we measured the 
temperature dependence of the specific heat. The clear anomaly at 17 K (figure 
[3]) and the linear behavior of C/T 2 (T) below this temperature are typical for 
antiferromagnets [38] . Thus, we interpret this as a transition to an AFM ordered state 
(Tn = 17 K). Remarkably, the anomaly does not shift nor decrease in amplitude in 
magnetic fields up to fioH = 9 T. Prior to the analysis of the magnetic behaviour 
above Tn, the specific heat should be decomposed into the magnetic contribution 
(which reflects the spectrum of magnetic excitations) and the phonon contribution (the 
spectrum of lattice vibrations). This decomposition is reliable only if the overlap of 
the two spectra (magnetic excitations and phonons) is relatively small (see [39] for an 
example). As the phonon contribution increases on temperature, the decomposition is 
possible for systems with weak magnetic couplings (J^- < 10 K). As we obtained from 
our susceptibility data, the energy scale of J in CuSe 2 5 is about 165 K. Thus, for a ID 
Heisenberg chain we expect the maximum of the magnetic specific heat at 0.48 J [36], i.e. 

§ A simultaneous fit of both curves implying the same J value in the whole temperature range down 
to the phase transition yields considerable deviations from experiment. The origin of this deviation is 
discussed below in the text. 



fitted both curves 




temperature window. As a result, 
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close to 80 K. At this temperature, the phonon contribution to the specific heat strongly 
dominates over the magnetic contribution. As a result, the experimental curve has no 
visible features in the vicinity of 80 K. For systems with large couplings (J^ > 10 K), 
the most accurate way to account for the phonon part is to measure an isostructural 
non-magnetic reference system (see jlO] for an example). In case of CuSe20s, it is not 
possible, as ZnSe20s has a different crystal structure [H], and thus a different phonon 
spectrum. Therefore, the specific heat data provide a clear evidence of an AFM ordering 
at 17 K but do not allow an independent justification of the temperature scale for the 
leading magnetic interactions. 



» n H=9T 
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Figure 3. C p /T 2 of CuSe20s as a function of temperature and magnetic field. The 
Neel temperature is marked with a dashed line. Inset (ordered phase region): the 
Cp^T 3 behaviour predicted by theory is complicated by a clear kink at 7 K. 



The last remark concerns a pronounced kink at 7 K (figure El inset), i.e. the 
region of the ordered phase. The kink is stable at least up to fi^H = 9 T and thus not 
related to defects. Intriguingly, a similar feature has been observed for a related system 
Bi2Cu04 (figure 3 in [12]) favouring the intrinsic nature of the kink rather than a sample 
dependent effect. For the magnetic contribution to the specific heat, such features have 
been proposed to mark a dying out of high frequency spin wave modes [3B]- To elucidate 
this unusual feature, further experimental studies on CuSe 2 5 and similar systems as 
well as a careful theoretical analysis should be carried out. 



4-2. Microscopic model 

We start in our microscopic analysis with band structure calculations performed in 
the local density approximation (LDA). LDA yields a valence band of about 9 eV width 
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formed mainly by Cu 3d, O 2p and Se 4p states (figurelH right panel). The well-separated 
double-peak at the Fermi level e-p contains two narrow, half-filled bands (figure HI left 
panel). The width of this antibonding band complex (0.85 eV) is in between the widths 
of the same complex in Bi 2 Cu0 4 (1.05 eV [31]) and Sr 2 Cu(P0 4 ) 2 (0.65 eV 0). The LDA 
yields a metallic GS, contrary to the experimentally observed insulating behavior. This 
discrepancy is caused by the underestimate of strong on-site Coulomb interactions of the 
Cu 3d electrons. Nevertheless, LDA reliably yields the relevant orbitals and dispersions. 
Thus, we have a closer look to the band complex at Ep. The two bands, relevant for 
the low-lying magnetic excitations, are related to the antibonding dp a orbital of a Cu0 4 
plaquette, i.e. the antibonding combination of Cu 3d x 2_ y 2 and O 2p a states (orbitals 
are denoted with respect to the local coordinate system). The antibonding dp a orbital 
is well separated (AE ~ 0.5 eV) from the lower lying Cu 3c? and O 2p states. Thus, 
the most efficient way to describe the electronic structure is to construct an effective 
one-band tight-binding (TB) model (one band per plaquette), parametrized by a set of 
electron transfer integrals tij. The correlation effects, insufficiently described by LDA 
and thus by the TB model, are accounted for by adopting a corresponding Hubbard 
model mapped subsequently onto a Heisenberg model (this mapping is valid for spin 
excitations in the strongly correlated limit at half-filling, both well justified for undoped 
cuprates with small magnetic exchange). 

Prior to numerical calculations, we compare the dispersions of the two well 
separated bands at (figure [5]) to dispersions of the corresponding antibonding dp a 
complexes of Sr 2 Cu(P04)2 (figure 2 in [8]) and Bi 2 Cu04 (figure 5 in (3lj). Here, a 
close similarity of CuSe 2 5 and Sr 2 Cu(P0 4 ) 2 is revealed: both band structures have 
a dominating dispersion along the chain direction (for CuSe 2 Os, this is the c-axis in 
figure [1] and T — Z region in figure H]) and a weaker dispersion in other directions, unlike 
Bi 2 Cu0 4 , where the dispersions along different directions in the fc-space are comparable, 
indicating a 3D behaviour. 

For a quantitative analysis, we constructed an effective one-band TB Hamiltonian 
and determined the set of transfer integrals in order to get the best least-squares 
fit to the two LDA bands crossing ep. As an alternative approach, we used a Wannier 
functions (WF) technique which implies the construction of WF for the Cu 3d x 2_ y 2 
antibonding state, relevant for the magnetism, and the calculation of the overlap of 
the WF. The results of the latter method are affected by the overlap of the relevant 
3d x 2_ y 2 antibonding state with other states. Therefore, for perfectly separated bands 
as in CuSe 2 05, both methods should yield the same results within numerical accuracy. 
The numerical evaluation supports this statement: the difference between the transfer 
integrals obtained by the WF method and by the TB fit is tiny and does not exceed 
2 meV for individual values (the mean value for all is 0.2 meV). This deviation can 
be considered as an error margin for the mapping procedure. Thus, for isolated bands 
the WF method should not be regarded more accurate than a direct TB fit, but rather 
as an independent alternative procedure pUl ELTJ 02] • The agreement of the results using 
the two independent mapping methods reflects the applicability of an effective one-band 
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approach. 

The resulting set of the transfer integrals (table [21 first column) yields perfect 
agreement with the LDA bands (figure [5]) IjjJ The hopping paths corresponding to the 
leading terms are shown in figure [5] (right panel) . 




Figure 4. LDA band structure (left panel) and density of states (right panel) of 
CuSe 2 5 . Notation of fc-points: r=(000), X=(f00), S=(f f 0), Y=(0f 0), Z=(00^), 
R=(^),A=(Sff),M=(0f^). 




r x 




y r z 
wave vector 



M 



Figure 5. Left panel: the tight-binding fit (circles) to the LDA band structure 
(antibonding dp a band, solid line) Right panel: the superexchange paths for the leading 
transfer integrals. The projection of the structure is the same as in Fig. [TJ 



We find that the leading couplings in CuSe2 0s are the NN intra-chain coupling 
t c = 165 meV and one of the short inter-chain couplings t a b ~ 45 meV. The corresponding 

|| To check the results for consistency, we have neglected all iy smaller than 10 meV and repeated the 
fitting. The difference of leading terms in both approaches did not exceed 10%. 



12 




Figure 6. Wannier functions for the Cu 3d x 2_ y 2 orbital. Colours represent the sign 
of a Wannier function. Left panel: the Cu 3d x 2_ y 2 Wannier function plotted on top of 
a CuC>4 plaquette, visualizing the antibonding combination of Cu 3d x 2_ y 2 and O 2p a 
states, relevant for the magnetism. Central panel: the overlap of two Wannier functions 
centered on the neighbouring Cu atoms (corresponds to the nearest neighbour intra- 
chain coupling t c ). Note that the neighbouring plaquettes are tilted which leads to 
a sizable 7T-overlap of the Wannier functions, and hence O 2p wave functions of the 
neighbouring plaquettes, allowing for a considerable ferromagnetic contribution to the 
magnetic exchange. Right panel: the overlap of the Wannier functions corresponding 
to the leading inter-chain coupling t a b- 

Wannier functions are pictured in figure El The value of the largest (NN intra- 
chain) coupling in CuSe20s is slightly larger than the corresponding coupling in 
Sr2Cu(PC>4)2 (135 meV [8]). The difference in the largest inter-chain term is more 
pronounced: the size of the inter-chain coupling in CuSe20s (45 meV) is considerably 
higher than in Sr 2 Cu(P04)2 (16 meV j8j). Even more important is the difference in 
the specific coupling geometry — whether it is constructive towards the long-range 
ordering or not. As we stated while comparing the crystal structures of CuSe 2 5 and 
Sr2Cu(PC>4)2 (see section [3]), in both systems there are two short inter-chain coupling 
paths. The corresponding couplings are identical (symmetry-related) in Sr 2 Cu(P04)2, 
but independent (and in fact, considerably different) in CuSe20s. The TB analysis 
reveals that only one (t a b) of the two NN inter-chain couplings is relevant for CuSe20s 
(table [2]). Consequently, the essential difference between the two systems can be best 
understood in terms of the spin lattices that are formed by the strongest intra-chain and 
inter-chain couplings, as depicted in figure [71 In Sr 2 Cu(P0 4 )2, three relevant couplings 
(the intra-chain NN coupling and two identical inter-chain couplings) are arranged on an 
anisotropic triangular lattice (figured left panel). By switching off one of the inter-chain 
couplings, the topology of the relevant couplings changes, and the system is described 
by two couplings forming an anisotropic square lattice (figure right panel). The 
main difference between the two topologies is that in Sr 2 Cu(P04)2 the competition of 
relevant couplings, which can not be simultaneously satisfied, leads to strong magnetic 
frustration, while in CuSe205 the inter-chain couplings are not frustrated. The lifting 
of frustration in CuSe20s has a remarkable influence on the physical properties as will 
be discussed below. 

The calculated transfer integrals provide valuable information on the coupling 
regime. To include the missing Coulomb interaction U e s, as described in the section El 
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Figure 7. Topology of inter-chain couplings in Sr2Cu(P04)2 (left panel) and CuSe20s 
(right panel). Red lines denote the AFM coupling. Bold red lines highlight the chains. 
In Sr2Cu(PC>4)2, the intra-chain and the two equivalent inter-chain couplings form 
an anisotropic triangular lattice. In CuSe20s, there is only one relevant inter-chain 
coupling resulting in an anisotropic square lattice geometry of intra-chain and inter- 
chain couplings. The former geometry leads to magnetic frustration, while the latter 
is not frustrated (see insets). 



we can use the TB model to construct a Hubbard model and map the latter onto 
a Heisenberg model to obtain the antiferromagnetic (AFM) exchange from J^ FM = 
4t| /t/ e if. Using the same representative U e g = 4.5 eV as for Sr 2 Cu(P04)2 [8], we obtain 
jAfm _ 235 x f or the NN intra-chain exchange and J^f M — 27 K for the largest inter- 
chain exchange. Other couplings yield values of AFM exchange less than 1.5 K (table El 
second column) and will be neglected in further discussion. 



Table 2. Leading transfer (first column) and exchange integrals (last column) of 
CuSe205. The AFM exchange (second column) is calculated via mapping the transfer 
integrals onto an extended Hubbard (U s — 4.5 eV) and subsequently onto a Heisenberg 
model. The total exchange is taken from LSDA+J7 total energy calculations of 
supcrcclls. The ferromagnetic exchange J^ M is evaluated as the difference between 



path %/meV J;f M /K j£ M /K J„/K 



x c 


166 


285 


-120 


165 


X a b 


51 


27 


-7 


20 


v2c 
A b 


11 


1.5 





1.5 


X ab 


10 


1 





1 


x c b 


10 


1 





<1 


v2c 
A a6 


7 


0.5 





<1 



The calculated leading magnetic exchange J A = 285 K is considerably larger 
than our estimate from the Bethe ansatz fit («155 K) based on experimental x{T) 
data. Moreover, it is larger than the corresponding exchange integral in Sr2Cu(P04)2 
(187 K). This discrepancy originates from ferromagnetic (FM) contributions to the total 
magnetic exchange, which are neglected in the mapping procedure. For the NN exchange 
J c , we expect a considerable FM contribution originating from the overlap of O 2p wave 
functions of neighbouring plaquettes. Due to a dihedral angle = 64° between the 
neighbouring plaquettes, this overlap has a sizable 7r contribution (this can be seen 
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in the WF in figure El central panel) leading to a Hund's rule (FM) coupling. For 
the leading inter-chain coupling, the FM contribution is expected to be small due to a 
predominantly a overlap of 2p wave functions (figure EJ right panel). 

To get a numerical estimate for the FM contribution, we perform total energy 
calculations for various spin patterns of magnetic supercells using the LSDA+f/ method. 
The method is rather sensitive to the U d value. As Ud = 6.5 eV yields agreement 
between the calculated and experimentally measured exchange integrals of the well 
studied La 2 Cu04 and CuGeC>3, we adopted this value in the calculations for CuSe 2 O5 
The supercell method has limitations set by numerical accuracy for the small exchange 
integrals (J^- < 1 K) and the size of the required supercells. In our case, we constructed 
supercells and spin patterns that yield all exchange couplings which were found to 
be relevant from the TB analysis. The resulting total energies are mapped onto a 
Heisenberg model, which is parameterized by the total exchange integrals (table EJ last 
column) containing both AFM and FM contributions. Thus, by subtracting the AFM 
part J(f M from the total exchange J^, the FM contribution J? M can be estimated 
(table [21 fourth column). 

In general, LSDA+t/ calculations yield a reliable estimate for exchange integrals [HI 
[T3J, EH E9] . This reliability holds for CuSe 2 5 : we obtain J c = 165 K and J ab = 20 K, in 
almost perfect agreement with the estimates from magnetic susceptibility. In accordance 
with our expectations, J™ = —120 K has a considerable contribution to the total 
exchange J c , while J™ = —7 K yields a smaller correction to the J a b value. We should 
note that the large ferromagnetic contribution J™ = —120 K may originate, in addition 
to the mentioned 7r-overlap of O 2p wave functions, also from a destructive interference 
of coupling paths [44] or by a strong coupling to ligands [45] . Which of these mechanisms 
plays a leading role in CuSe 2 5 is an open question. This issue is, however, beyond the 
scope of the present paper and needs further theoretical investigation. 

Though the calculated J c value (165 K) is very close to the estimate from the Bethe 
ansatz (157 K), we decided to check the exchange integrals for consistency by performing 
additional calculations for Ud = 6.0 eV and Ud = 7.0 eV. Besides the expected change 
of exchange integrals (0.5 eV increase of the Ud results in about 20% decrease of 
and vice versa), we found that the ratio a = J a b/Jc of the leading exchange integrals 
(a = 0.121 for U d = 6.0 eV, a = 0.129 for U d = 6.5 eV and a = 0.136 for U d = 7.0 eV) 
is rather stable with respect to the Ud value. 

Thus, the consideration of the FM contribution yielded a valuable improvement of 
the energy scale comparing to the AFM exchange values J^ FM and J^ FM , but the ratio 
a of the two couplings, that is the most relevant for the magnetic ground state, stays 
almost unchanged. Moreover, this ratio is stable with respect to the model parameters 
U c ff and Ud, leading to a very reliable physical picture: CuSe 2 5 can be described as 

% It is worth to note, that the Ud parameter is not universal and depends on a calculational scheme 
and consequently on the basis set implemented in a code. Thus, different Ud values adopted in this 
work for CuSe 2 5 (U d = 6.5 eV, the code fplo version 7.00-27) and for Sr 2 Cu(P0 4 ) 2 (U d = 8.0 eV [8], 
the code fplo version 5.00-18) originate from the different basis used in the codes. 
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a quasi- ID system with AFM chains characterized by an NN intra-chain exchange of 
165 K. Each chain is coupled to two neighbouring chains by the non-frustrated inter- 
chain exchange of 20 K (one order of magnitude smaller that the intra-chain coupling). 

4-3. Simulations 

Turning back to the discussion about the ordering temperature (section [I]), it is 
reasonable to point out the advantages of CuSe20s as a model system. First, we have 
evidence from both theory and experiment that the system is mainly ID. Secondly, the 
microscopic analysis revealed that the NNN intra-chain coupling is practically absent 
leading to a valuable simplification for a theoretical analysis. Finally, there is only 
one relevant inter-chain coupling. The fact that this latter coupling is not frustrated 
allows to use the powerful quantum Monte Carlo (QMC) method for a simulation of 
thermodynamical data with a subsequent comparison to the experimentally measured 
curves. The results of the simulations are given in figure M in comparison with the 
Bethe ansatz fits (where the inter-chain coupling is neglected) . Obviously, the inclusion 
of the inter-chain coupling yields only a tiny improvement with respect to the Bethe 
ansatz fits. This fact demonstrates a posteriori the importance of a microscopic model 
for systems like CuSe20s: apart from the microscopic modelling, there is no reliable 
way to account for the small inter-chain coupling directly from measurements of the 
paramagnetic susceptibility. 




100 200 300 



T(K) 



Figure 8. Comparison of Bethe ansatz and quantum Monte Carlo (QMC) fits to the 
experimental magnetic susceptibility. The temperature-independent contribution \o 
in both fits was set to zero. 
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4-4- Estimation of the Neel temperature 

To benefit from the unique combination of the simple microscopic picture and the 
experimentally well-determined AFM ordering temperature Tn, we make an attempt to 
estimate Tn from an available simplified theory using the calculated exchange integrals. 
Still, there are three problems, intrinsic for quasi-lD systems, to be accounted for. The 
first is the spatial anisotropy of exchange couplings present in a system. This problem 
is resolved in CuSe20s only partially. On one hand, there is only one relevant inter- 
chain coupling for every pair of neighbouring chains, but on the other hand, it couples a 
certain chain with only two of four neighbouring chains. Thus, the couplings to the other 
two chains are considerably smaller, resulting in the (spatial) exchange anisotropy. The 
second problem is the anisotropy in the spin space. In our microscopic approach, we used 
the isotropic Heisenberg model where this anisotropy is neglected. Although we observe 
a remarkable agreement between the microscopic model and the macroscopic behaviour, 
the spin anisotropy is present, as evidenced for instance by the strong dependence of 
the g- factor on the orientation of a magnetizing field (figure [2]). In a common sense 
approach, the (spatial) exchange anisotropy is expected to lower T N , while the spin 
anisotropy raises it. In each system, these two effects are balanced. Attempts to find 
a suitable description for this balance were made in a number of advanced theoretical 
studies based on a mean-field formalism [181 US], [21] [2l]- Still, the problem seems 
not to be resolved, since a considerable disagreement remains between the numerical 
results yielded by different theories (none of which is generally accepted) and, even more 
important, due to the third problem — the problem of inter-chain magnetic frustration, 
which has not been addressed, so far. 

Here, we make an empirical attempt to estimate how the frustration influences the 
magnetic ordering. For that purpose, we compare several quasi-lD magnetic compounds 
in a systematic way. The two well studied quasi-lD cuprates Sr 2 Cu03 and Ca2Cu03 
are commonly referred as model systems in most theoretical studies regarding the Tn 
problem. We should note that these two systems are essentially different from CuSe205 
due to the presence of corner-sharing chains of CUO4 plaquettes, which results in 
one order of magnitude larger NN coupling. In addition, the NNN coupling is not 
negligible [6]. Nevertheless, they are referred here for the sake of completeness. In 
both, Sr 2 Cu03 and Ca2Cu03, one of the relevant inter-chain couplings is frustrated. To 
calculate Tn, we use formulas given by Schulz [19] . Irkhin and Katanin [21] and Yasuda et 
al. [23]. The calculated Tn are given in the three last columns of (table [3]). Disregarding 
the method used, the calculated Tn considerably overestimate the experimental values 
for Sr 2 Cu03 and Ca2Cu03 (table El fifth column). 

A theoretical approach is expected to work better for systems with more pronounced 
ID nature — e.g. Sr2Cu(P04)2 and K2CUP2O7. The structural peculiarities of these 
systems were discussed in section |3j The main issue here is the frustration caused 
by the leading inter-chain coupling. For Sr 2 Cu(P04)2, theory predicts an ordering 
temperature Tn two orders of magnitude larger than the experimentally observed value. 
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Table 3. Exchange integrals (columns 2 and 3) together with experimental (column 
5) and theoretically calculated (columns 7-10) ordering temperatures Tn for quasi- 
onc-dimensional cuprates (3D magnet Bi 2 Cu04 has a similar structural motive and 
was added for completeness). The ordering temperatures Tn were calculated using 
formulas from references [19l [21], [23] and the exchange integrals from columns 2 (nearest 
neighbour intra-chain coupling Ji) and 3 (leading inter-chain coupling Jj_). Note that 
the first four systems are frustrated due to inter-chain couplings, while in the last two 
the inter-chain couplings are not frustrated. 



compound 


Ji/K 




Ref. 


Tn exp./K 


Rcf. 


m 


Tn 

m 


calc./K 

[23] [31] 


Sr 2 Cu0 3 


2200 


9 


i 


5 


m 


28 


22 


21 


Ca 2 Cu0 3 


1850 


42 




9 


m 


115 


91 


85 


Sr 2 Cu(P0 4 ) 2 


187 


3 


i 


0.085 


m 


8.5 


6.7 


6.3 


K 2 CuP 2 7 


196 


0.25 


m 


< 2 


m 


0.9 


0.7 


0.6 


CuSe 2 5 


165 


20 




17 




23 


18 


17 


Bi 2 Cu0 4 


10 


G 


m 


42 


[32] 






47 



(Unfortunately, experimental low-temperature data are not available for K2CUP2O7.) 
This huge discrepancy is in sharp contrast with the situation for the 3D magnet Bi 2 Cu04, 
for which the theoretical estimate coincides with the experimental value within the error 
bars (table El last row) El 

In CuSe205, each chain is strongly coupled only with two of four neighbouring 
chains (unlike Sr 2 Cu03, Ca 2 Cu03, Sr 2 Cu(P04)2 and K2CUP2O7 with coupling to four 
neighbouring chains). There is no unique way to take this feature into account. A simple 
approximation is to take the arithmetic average, which yields an effective inter-chain 
coupling value J± = J a b/2^ Using this value, theory yields a perfect agreement with 
experimental value (tabled fifth row). 

Obviously, the existing models describe the magnetic ordering in CuSe20s much 
better than in Sr 2 Cu0 3 and Ca 2 Cu0 3 , and especially in Sr 2 Cu(P0 4 )2. Despite our crude 
way of accounting for spatial exchange anisotropy and the neglect of spin anisotropy, 
the theoretical estimate of T N for CuSe 2 5 is in surprisingly good agreement with 
the experimental value. Though in general the magnetic ordering is affected by the 
spin anisotropy, CuSe20s yields empirical evidence that for systems with a small spin 
anisotropy the isotropic model provides a rather accurate estimate of T N . Thus, it 
is unlikely that the disagreement between theoretical and experimentally observed Tn 
values for Sr 2 Cu(P04)2 originates from the neglect of spin anisotropy effects. 

Finally, only the magnetic frustration is left to be a possible reason for a huge 
discrepancy between theory and experiment. Our analysis reveals that frustrated inter- 
chain couplings play a crucial role for the magnetic ordering. This fact explains why 

+ For the theoretical estimation of Tn, the formula 7 from [31] was used. 

* Alternatively, the geometrical averaging can be used. This approach yields a correct limit with 
respect to the Mermin- Wagner theorem (zero ordering temperature for ID and 2D systems). Then, 
calculational schemes from [19], [21] and [23] yield Tn values of 14 K, 11 K and 10 K, respectively. 
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theoretical schemes fail to predict T N for frustrated systems. 

To illustrate the influence of frustration, we use a simple formula from the spin 
wave theory in a random phase approximation, which connects Neel temperatures 
for two compounds A and B with the values of exchange integrals: T^/T^pa 



y Ji-J^l 'y Jf 'J® [6]. Using the values of exchange integrals for Sr 2 Cu(P0 4 )2 and 

CuSe205 and the experimental Neel temperature for Sr 2 Cu(P04)2 (table [3]), we obtain 
Tn pa 0.146 K for CuSe205, almost 120 times smaller than the experimental value. 

In the existing theoretical approaches, a parameter controlling the frustration 
caused by inter-chain couplings is missing. Therefore, new theories which would treat 
magnetic frustration as one of the key issues for the magnetic ordering, are needed. On 
the other hand, there is a lack of information from the experimental side, resulting in 
a very limited number of systems that challenge the theoretical predictions. Besides 
CuSe205, an almost perfect model system, synthesis and investigation of new systems 
with similar crystal chemistry are highly desirable. 

5. Summary and outlook 

The class of quasi- ID magnets attracts much attention as a field of search for prominent 
models and a playground for modern theories. Recently, by studying the magnetic 
properties of Cu 2+ phosphates, several systems of this class were found to exhibit 
the physics of a Heisenberg chain model. In these materials, the remarkable one- 
dimensionality and the absence of long-range intra-chain interactions are ruled by a 
unique arrangement of magnetically active CUO4 plaquettes: they form edge-sharing 
chains where every second plaquette is cut out. The chains are well separated by 
alkaline or alkaline earth cations (K, Sr). The magnetic susceptibility of these systems is 
perfectly described by the Bethe ansatz, which provides an exact solution for the nearest 
neighbour (NN) antiferromagnetic spin-1/2 Heisenberg chain. At the same time, the 
ordering temperature Tn of the systems reveals a fundamental disagreement between 
theory and experiment. Unfortunately, the range of available experimental studies of 
these systems is rather limited, as the materials are presently available only as powders. 

Therefore, we have synthesized CuSe20s — a system implying a similar, isolated 
arrangement of neighbouring Cu0 4 plaquettes (but tilted with respect to each 
other, unlike Cu 2+ phosphates) and allowing for a growth of high quality single 
crystals. Thermodynamic measurements reveal a quasi- ID behaviour with a leading 
antiferromagnetic coupling of about 160 K (obtained from the Bethe ansatz fit for 
the magnetic susceptibility). The system orders antiferromagnetically at 17 K, as 
evidenced by magnetic susceptibility and specific heat data. A microscopic analysis 
based on the results of DFT calculations reveals that CuSe20s can be described in good 
approximation by only two relevant exchange integrals: NN intra-chain (J c = 165 K) 
and the leading inter-chain coupling (J a b = 20 K). The theoretical estimate of 
the ordering temperature Tn is in perfect agreement with experimental value. This 
remarkable agreement is in sharp contrast with a huge overestimate of Tn for Cu 2+ 
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phosphates yielded by a formal application of the same theory. To reveal the origin of 
this difference on empirical grounds we analyzed systematically the factors affecting T^. 
Beyond the influence of the spatial exchange anisotropy and the spin anisotropy, we 
emphasize the role of the magnetic frustration due to equivalent inter-chain interactions 
in the latter compounds. Comparing theoretical and experimental data for related 
systems, we show that inter-chain frustrations have a crucial influence on Tn and are 
likely the main cause for the failure of any theory which ignores them. 

For an outlook, we propose further experimental studies (for instance, Raman 
spectroscopy and inelastic neutron scattering) in order to benefit from the availability 
of single crystals of CuSe20s. Especially, additional experimental data are required in 
order to understand the nature of the specific heat anomaly at 7 K. Secondly, we hope 
that our work will inspire a directed search for new quasi-lD model systems. Last but 
not the least, we want to stimulate the development of more sophisticated theories for 
the estimation of T^. In particular, such theories should explicitly take into account the 
magnetic frustration arising from complex inter-chain interactions. The well understood 
system CuSe 2 5 could give valuable support for these theories. 

Acknowledgements 

We are grateful to P Scheppan for her help with the EDXS analysis and to R Cardoso- 
Gil and J A Mydosh for critical reading of the manuscript and valuable suggestions. The 
investigation was supported by the GIF (1-811-237.14103) and by the Emmy Noether 
program of the DFG. 

References 

[1] Hasc M, Tcrasaki I and Uchinokura K 1993 Phys. Rev. Lett. 70 3651-3654 

[2] Enderle M, Mukherjee C, Fak B, Kremer R K, Broto J M, Rosner H, Drechsler S L, Richter J, 

Malek J, Prokofieva A et al 2005 Europhys. Lett. 70 237-243 
[3] Drechsler S L, Volkova O, Vasiliev A N, Tristan N, Richter J, Schmitt M, Rosner H, Malek J, 

Klingeler R, Zvyagin A A et al 2007 Phys. Rev. Lett. 98 077202 |arXiv:cond-mat/0701741| 
[4] Schmitt M, Malek J, Drechsler S L and Rosner H To be published 
[5] Bethe H 1931 Z. Phys. 71 205-226 

[6] Rosner H, Eschrig H, Hayn R, Drechsler S L and Malek J 1997 Phys. Rev. B 56 3402 

|arXiv:cond-mat/9704210| 
[7] Ami T, Crawford M K, Harlow R L, Wang Z R, Johnston D C, Huang Q and Erwin R W 1995 

Phys. Rev. B 51 5994-6001 
[8] Johannes M D, Richter J, Drechsler S L and Rosner H 2006 Phys. Rev. B 74 174435 

arXiV:cond-mat/0609430| 
[9] Belik A A, Azuma M and Takano M 2004 J. Solid State Chem. 177 883-888 
[10] Salunke S S, Ahsan M A H, Nath R, Mahajan A V and Dasgupta I 2007 Phys. Rev. B 76 085104 
[11] Rosner H, Schmitt M, Kasinathan D, Ormeci A, Richter J, Drechsler S L and Johannes M D 2009 

Phys. Rev. B 79 127101 

[12] Salunke S, Ahsan M A H, Nath R, Mahajan A V and Dasgupta I 2009 Phys. Rev. B 79 127102 
[13] Nath R, Kasinathan D, Rosner H, Baenitz M and Geibel C 2008 Phys. Rev. B 77 134451 
larXrV:0804.1262l 



20 



Farnell D J J and Parkinson J B 1994 J. Phys: Condens. Matter 6 5521-5532 
Okamoto K and Nomura K 1992 Phys. Lett. A 169 433-437 

Chitra R, Pati S, Krishnamurthy H R. Sen D and Ramascsha S 1995 Phys. Rev. B 52 6581-6587 
Tonegawa T and Harada I 1989 J. Phys. Soc. Japan 58 2902-2915 
Scalapino D J, Imry Y and Pincus P 1975 Phys. Rev. B 11 2042-2048 
Schulz H J 1996 Phys. Rev. Lett. 77 2790-2793 |arXiv:cond-mat /9604144] 

Essler F H L, Tsvclik A M and Delfino G 1997 Phys. Rev. B 56 11001-11013 



arXiv:cond-mat/9705196 



Irkhin V Y and Katanin A A 2000 Phys. Rev. B 61 675 7-6764|arXiv:cond-mat/9909257 
Bocquet M 2002 Phys. Rev. B 65 184415 |arXiv:cond-mat/0110429| 

Yasuda C, Todo S, Hukushima K, Alet F, Keller M, Troyer M and Takayama H 2005 



Phys. Rev. Lett. 94 217201 arXiv:cond-mat/0312392 



Todo S and Shibasaki A 2008 Phys. Rev. B 78 224411 [arXlE:0805.3097| 
Zvyagin A A and Drechsler S L 2008 Phys. Rev. B 78 014429 larX"iv:0806.1074l 
Kahn O, Verdaguer M, Girerd J J, Galy J and Maury F 1980 Solid State Commun. 34 971-975 
Koepernik K and Eschrig H 1999 Phys. Rev. B 59 1743-1757 
Pcrdcw J P and Wang Y 1992 Phys. Rev. B 45 13244-13249 
Eschrig H, Koepernik K and Chaplygin I 2003 J. Solid State Chem. 176 482-495 
Albuquerque A F, Alet F, Corboz P, Dayal P, Feiguin A, Fuchs S, Gamper L, Gull E, Giirtler S, 

Honccker A et al 2007 J. Magn. Magn. Mater. 310 1187-1193 
Janson O, Kuzian R O, Drechsler S L and Rosncr H 2007 Phys. Rev. B 76 115119 
Petrakovskii G A, Sablina K A, Pankrats A I, Vorotinov V M, Furrer A, Roessli B and Fischer P 

1995 J. Magn. Magn. Mater. 140-144 1991-1992 
Meunier G, Svenssen C and Carpy A 1976 Acta Crystallogr. B32 2664 
Becker R and Berger H 2006 Acta Crystallogr. E62 1256 1257 

Janson O, Richter J and Rosncr H 2008 Phys. Rev. Lett. 101 106403 larXiv:0806.1592l 
Johnston D C, Kremer R K, Troyer M, Wang X, Kliimper A, Bud'ko S L, Panchula A F and 

Canfield P C 2000 Phys. Rev. B 61 9558-9606 (see table I, fit 2) 
Goodcnough J B 1963 Magnetism and the Chemical Bond (New York: Interscience- Wiley) 
St0len S, Grande T and Allan N L 2004 Chemical Thermodynamics of Materials (Chichester: John 
Wiley & Sons, Ltd) 

Nath R, Tsirlin A A, Rosner H and Geibel C 2008 Phys. Rev. B 78 064422 [arXlv:0803.3535l 
Kondo S, Johnston D C, Swenson C A, Borsa F, Mahajan A V, Miller L L, Gu T, Goldman A I, 

Maple M B, Gajewski D A et al 1997 Phys. Rev. Lett. 78 3729-3732 
Meunier G and Bertaud M 1974 Acta Crystallogr. B30 2840-2843 

Baran M, Gaidukov Y P, Danilova N P, Inushkin A V, Jedrzejczak A, Koksharov Y A, Nikiforov 
V N, Revcolevschi A, Szymczak R and Szymczak H 1999 J. Magn. Magn. Mater. 196-197 
532-533 

Tari A 2003 The specific heat of matter at low temperatures (London: Imperial College Press) 
p 148 

Feldkemper S, Weber W, Schulenburg J and Richter J 1995 Phys. Rev. B 52 313-323 
Geertsma W and Khomskii D 1996 Phys. Rev. B 54 3011-3014 

Kojima K M, Fudamoto Y, Larkin M, Luke G M, Merrin J, Nachumi B, Ucmura Y J, Motoyama 
N, Eisaki H, Uchida S et al 1997 Phys. Rev. Lett. 78 1787-1790 
[47] Bclik A A, Uji S, Terashima T and Takayama-Muromachi E 2005 J. Solid State Chem. 178 3461- 
3463 



